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The ability to represent intracellular biochemical dynamics via deterministic and stochastic 
modelling is one of the crucial components to move biological sciences in the observe-predict-control- 
design knowledge ladder. Compared to the engineering or physics problems, dynamical models in 
quantitative biology typically dependent on a relatively large number of parameters. Therefore, the 
relationship between model parameters and dynamics is often prohibitively difficult to determine. 
We developed a method to depict the input-output relationship for multi-parametric stochastic 
and deterministic models via information-theoretic quantification of similarity between model 
parameters and modules. Identification of most information-theoretically orthogonal biological 
components provided mathematical language to precisely communicate and visualise compensation 
like phenomena such as biological robustness, sloppiness and statistical non-identifiability. A 
comprehensive analysis of the multi-parameter NF-^B signalling pathway demonstrates that the 
information-theoretic similarity reflects a topological structure of the network. Examination of the 
currently available experimental data on this system reveals the number of identifiable parameters 
and suggests informative experimental protocols. 

Supplementary Information available at: 
http://www. ippt .gov.pl/~mkomor/redundancySI .pdf 



Last decades accumulated sufficient evidence that a num- 
ber of biological phenomena, in particular those related 
to intra-cellular dynamics, noise management, biochem- 
ical signalling cannot be understood by intuition alone 
and require mathematical formalism to explain and sum- 
marise available data. Expectably, mathematical mod- 
elling will help in prediction, control and design of bio- 
chemical networks. Therefore adaption of conventional 
modelling techniques is required to suit the specificity of 
these problems. Models of biochemical dynamics are dif- 
ferent from classical models of engineering and physics in 
a number of ways. Primarily they involve substantially 
larger relative number of parameters compared to avail- 
able data size. This challenge has given rise to a number 
of approaches aimed at improving our ability to develop, 
verify and apply mult i- parameter mechanistic models of 
such systems. We can loosely group these methods into 
those aimed at determining model sensitivities to param- 
eter values pj[3] , tools to estimate rate parameters ]4][8j , 
and techniques focused on maximisation of the informa- 
tion content of the experimental data [9HT2]. The in- 
put (parameters) - output (dynamics) dependencies is 
the main considered object of the above methods. The 
concept of information, which in the Fisher sense is a 
sensitivity of an output to parameters, establishes a nat- 
ural language to communicate a number input-output 
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phenomena. Sensitive parameters exert strong impact 
on output and therefore are relatively easy to infer. In 
consequence, when aiming at parameters estimation, ex- 
perimental settings, which render model parameters sen- 
sitive, should be searched. A number of studies have re- 
ported the intrinsic feature of dynamic multi-parameter 
models of biochemical dynamics to be sensitive only to a 
small number of linear combinations of parameters [2| [T3l- 
H5] . The developed methodology substantial enriched our 
repertoire of techniques to investigate input-output rela- 
tionship in multi-parameter models [THSl [TTJ HMST] . In 
this paper we built upon these findings to take a com- 
prehensive view at the problem of sensitivities in multi- 
parameter models. A notion of functional redundancy 
between individual parameters and their groups (mod- 
ules) is introduced with Shannon Information being a 
measure of its strength. As a result we propose a natural 
and general mathematical language to precisely commu- 
nicate and visualise all types of compensation like phe- 
nomena i.e. multi-parameter sensitivity, biological ro- 
bustness, sloppiness and statistical non-identifiability. It 
allows for a more insightful interpretation of sensitivity 
coefficients, detection and elimination of non-identifiable 
parameters and guided design of experiments aiming at 
maximising the number of identifiable parameters. We 
also find two efficient ways to evaluate functional redun- 
dancy. One is based on the Fisher Information (FI), 
therefore is local in the parameters space and requires 
parameter values as input; second is local in the space of 
experimental results and is based on posterior distribu- 
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tion sampling. We also integrate the introduced redun- 
dancy measure with a hierarchical clustering algorithm 
to informatively represent redundancy structure in form 
of dendrograms so that functionally related / orthogonal 
biological components can be easily identified and con- 
clusions about sensitivity, identifiability, robustness can 
be made. 

The method and its underlying principles are very gen- 
eral and applicable to deterministic and stochastic mod- 
els. The limiting factor is the computational power. 
In this paper we focus mainly on the computationally 
least demanding scenario, deterministic model with avail- 
able parameter guesses, and report unprecedented in- 
sight about parameters redundancies and their conse- 
quences for systems biology modelling. The potential of 
our framework is demonstrated using an example of the 
NF-ftB signalling pathway. Interestingly, we find that re- 
dundant parameters describe modules that are close in 
the network topology. We analyse experimental proto- 
cols published in the literature [22-30 and find that 26 
out of 39 model parameters can be estimated from litera- 
ture data. In addition, we use our method to propose 10 
stimulation experimental protocols, which are expected 
to provide estimates of 7 parameters non-identifiable so 
far. In the Supplementary Information (SI) we also cal- 
culate redundancies for a JAK-STAT model [91I3T]. We 
consider a stochastic model and a deterministic model 
with unknown parameter values. 



I. FUNCTIONAL REDUNDANCY (METHODS) 

A typical biochemical kinetics model describes how 
the abundances of a set of k molecular entities, y = 
(?/i, yk), change with time, t. Most generally the dis- 
tribution over times, (ti,...,t n ), of abundances, Y = 
(2/(ti), y(t n )), can be written as a probability distribu- 
tion P(Y\0), where = is a vector containing 
model parameters. Often y(t) is modelled as a solution 
of ordinary or stochastic differential equation 

^=F{y,0)+t{t), (1) 

where £(t) constitutes a random perturbation in the 
stochastic setting and is not included = 0) in the 

deterministic regime. Practically, only certain compo- 
nents of y can be measured experimentally with certain 
measurement error e. Therefore we consider variable JT, 
which contains selected elements of Y and measurement 
noise. In the Bayesian setting, the distribution of data 
given parameters, P(X\0), together with the prior dis- 
tribution, P(0), define, through averaging over possible 
parameter values, the distribution of possible measure- 
ments, P(X). The uncertainty in the possible measure- 
ments can be formally quantified in terms of the Shannon 
entropy, H(X). The average reduction in entropy of X 
gained by knowing is given by the mutual information 



I(X,0) = H(X) - H(X\6). The entropy of X there- 
fore can be seen as resulting from uncertainty in and 
other sources: H(X) = I(X,0) + H(X\0). Assume de- 
composition of parameters vector = (0a, 0b) and set a 
component- wise independent prior P(0) = P(0a)P(0b)- 
The entropy H(X) can now be divided [20] into con- 
stituents resulting from components of 

H(x) = i(x,e A ) + i(x,e B ) + i(e A ,e B \x) (2) 

+ H(X\{6 A ,6 B }). 

The reduction of entropy resulting from the sole knowl- 
edge of either 9a or Ob is described by I (0a) and I (9b) 
respectively. Mutual information I(0a,0b\X) measures 
the part of entropy corresponding to the concurrent 
knowledge of 0a and 0b- Intuitively, if 0a and 9b have 
redundant role knowing one parameter will not reduce 
uncertainty about X as it will be reproduced by uncer- 
tainty in the remaining parameter. Therefore we propose 
to use I(0a,0b \X) as a natural measure of redundancy 
between model parameters. The redundancy has also a 
dual interpretation. The reduction in uncertainty of 9 
gained by observing data X decreases the entropies of 
9a and ^ to a degree dependent on their redundancy 
level measured by I(0a,0b\X). See Fig. [i] for illustra- 
tion and SI for technical details. 

Evaluation of I(0a,0b\X) in a general setting it is com- 
putationally prohibitively expensive for most models of 
realistic size (see SI for details). We show, however, how 
to evaluate I(0a,0b\X) efficiently in two relevant sce- 
narios. In the first one, we assume that guesses about 
true parameter values are available, in the second that 
experimental data can be used to generate posterior dis- 
tribution. Denote available parameter guesses as 0* . In 
this case the posterior distribution can be approximated 
using the bayesian asymptotic theory 

P{6\6*) oc exp(-^(0 - 9*)FI(9*)(9 - 8*f), (3) 

where FI is the Fisher Information matrix of the model 
P(X\0*) and N is the number of replicates of X consid- 
ered in the above posterior. If FI has the full rank then 
the above density describes the multivariate normal dis- 
tribution (MVN). In this case we found a straightforward 
and numerically stable method to evaluation mutual in- 
formation 

1 m 

I(9 A ,9 B \9*) = --Y J ^-P]), (4) 

i=i 

where pj are canonical correlations extracted from the 
FI(9*) and m is a minimum of lengths of 0a and 0b 
(see SI for derivation). Canonical correlations are de- 
fined as the maximal correlation between mutually or- 
thogonal linear combinations of parameters of 0a and 0b 
(see SI). Therefore we can easily quantify functional re- 
dundancy if FI can be computed. This can be efficiently 
done, virtually without any computational limitations, 
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FIG. 1. Illustrative description of the functional redundancy. The uncertainty about X arises from an uncertainty about 
parameters 6 — (6a, 6b) via model P(X\6) and distributions P(0a) and P(6b)- The part of the entropy of X that can be 
explained by 6 is the mutual information I(X,6), which arrises from its three constituents: part of the entropy which can be 
explained solely by 6a i.e. I(6a,X); solely by 6b i.e. I(6b,X); and only by concurrent knowledge of 6a, 6b i.e. I(6a,6b\X). If 
6a and 6b are highly redundant knowing one parameter at a time will not reduce uncertainty about X as it will be reproduced 
to by the uncertainty about the remaining parameter. Dually, the reduction in uncertainty about 6 resulting from data, X, 
does not directly translates into reduction of individual entropies of 6a or 6b if these are functionally redundant. 



for deterministic models and for models of moderate size 
in the stochastic setting [3]. The second computation- 
ally tractable case is the one with available experimental 
data, X*. Here mutual information I(6a,0b\X = X*) 
can be estimated based on a sample from a posterior 
distribution P(9\X*). Sampling can be performed using 
one of the available Monte Carlo approaches (see SI for 
details). We show that despite limitations to evaluate 
computational redundancy in a general setting an un- 
precedented insight can be gained by computing redun- 
dancy in the two above scenarios. For simplicity, from 
now on we write 1 (6a, 6b) instead of I(0a,0b\Q*) an d 
I(0 A ,0 B \X*). 

Clustering to reveal functional parameters redundancy. 
Being able to evaluate 1(6 a, fe) we can represent redun- 
dancy structure for a given model using statistical clus- 
tering tools. Among a number of clustering approaches 
[32j [33] we chose the hierarchical clustering, which is 
more appropriate to represent correlation structure with 
assignment to clusters having a secondary importance 
(see SI for details). In order to built an intuition how a 
model can be analysed using this approach we start with 
the toy model of gene expression. 

Gene expression. Assume that the gene expression pro- 
cess begins with the production of mRNA molecules (r) 
at rate k r . Each mRNA molecule may be independently 
translated into protein molecules (p) at rate k p . Both 
mRNA and protein molecules are degraded at rates j r 
and 7 P , respectively. We consider steady state behaviour 
of deterministic model. Therefore, we have the state vec- 
tor (r,p) = TfT^- This formula suggests that pa- 



rameters k r and 7 r are entirely redundant as is the pair 
k p and j p . Parameters pairs (fc r ,7r) an d (k p ,j p ) are 
only partially redundant as the ratio k r j^ r has identi- 
cal impact on protein level as the ratio k p /j p and the 
latter does not impact mRNA level. This redundancy, 
reproduced by clustering algorithm, is visualised by the 
dendrogram in Figure [2j Parameters k r ,j r and k p ,j p 
are first linked reciprocally at hight and both pairs 
are linked together at non-zero high. We plot the link- 
ages at height — ^ Ej=i — pf)i (where m is a size 
of a smaller cluster) compared to the maximum of sons 
height. Canonical correlation gives a clear interpreta- 
tions to linkages heights, which is between and 1. The 
heatmap (top left corner) is the normalised FI matrix, 
I FI ij(0 ) ■ g ens jtivity coefficients, FIu(Q), are 

plotted in the bottom left corner. The steady state be- 
haviour of this simple system is robust to mutually com- 
pensating perturbations of production and degradation 
rates. The compensation renders two parameters identi- 
fiable, one of each of the pairs (k r ,j r ) and (k p: j p ). This 
redundancies can be removed by manipulating the initial 
condition and observing system temporal dynamics. 

The p53 system. In the above simple linear model re- 
lationship between behaviour and parameters can be un- 
derstood well without our information theoretic analysis. 
In a more complex, though still simple example, model 
of the p53 signalling system, a feedback loop between the 
tumor suppressor p53 (x p ), the oncogene Mdm2 mRNA 
(xq) and Mdm2 (xi), generates oscillation in response to 
DNA damage. Hence, the relations between parameters 
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FIG. 2. Functional redundancy for the simple model of gene 
expression. (A) Dendrogram. Linkages are plotted at height 
— m Sj=i l°9(X ~ Pj)- (B) Sensitivity coefficients: diagonal 
elements of the FI matrix. (C) Normalised FI matrix. (D) 
Sensitivity spectrum i.e log-eigenvalues of the FI matrix. Pa- 
rameters used: k r = 100 , k p = 2, 7 r = 1.2, 7 P = 1. 



is rather obscure. The deterministic version of the model 
is formulated as follows [34] 

xo = PyX p - a x 
x\ = ao^o — OL y x\. 



QiJzXy 



(5) 



Our approach can provide a comprehensive qualitative 
summary of parameters redundancy given that the pa- 
rameter values are provided. We use values published 
in [34]. The dendrogram in Fig. [3]A. predicts that of all 
parameters 6 = (f3 X: a X: a^, fc, (3 y , ctQ, a y ) the ao,a x ex- 
ert most similar impact on model dynamics and b x ,a y 
are 'least similar' ones. This is confirmed by plotting 
derivatives of model trajectories (Fig. |3^) with respect 
to these parameters. Conclusions regarding model ro- 
bustness and identifiability can be easily analogised to 
the previous example. 
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FIG. 3. Redundancy analysis of the p53 model. (A) Den- 
drogram. (B) Sensitivity coefficients. (C) Normalised FI 
matrix. (D) Log - eigenvalues of the FI matrix. (E) Deriva- 
tives of the model trajectories with respect to selected pa- 
rameters. Derivatives of ao,a x ,ak exhibit almost identical 
pattern, which is distinct from of this of a y , b x . Parameters 
used: /3 X 90, a x .002, a k 1.7, k 0.01, f3 y 1.1, a 



0.8, a v 



0.8. 



II. RESULTS 

The NF-kB system. Here we demonstrate the potential 
of our method using a deterministic model of the NF-^B 
signalling pathway, one of the key components control- 
ling innate immune response. In particular we address 
four important problems: 1) identify most information- 
theoretically orthogonal components of the network that 
control its behaviour; 2) characterise robustness proper- 
ties of the pathway; 3) analyse all published experimen- 
tal protocols to asses which parameters can be estimated 
from the available data; 4) propose stimulation proto- 
cols, which can increase the number of identifiable kinetic 



constants. The model considered here, proposed in [13] 
and further developed in [30] represents activation of NF- 
ftB induced genes in response to stimulation by TNF-a 
(the pro-inflammatory cytokine). It involves 39 parame- 
ters and 19 variables and encapsulates typical features of 
systems biology models that challenge current modelling 
techniques. 

Redundant control of the response to TNF-a stimulation. 
Visualisation of the redundancy structure enables identi- 
fication of network components that contain most redun- 
dant parameters and those most orthogonal that con- 
trol network dynamics. We first analysed scenario where 
behaviour of the system, X, is defined by trajectories 
of all model variables normalised by their maxima mea- 
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sured with normal measurement error e ~ MVN(0, 1) 
(Xi = Yi/max(Yi) + e). We considered reminiscence 
of physiological TNF-a stimulation profile and assumed 
TNF-a concentration to increase and drop after an inter- 
mediate plateau (see Figure IB of SI). The constructed 
redundancy dendrogram (Fig. [4| interestingly indicates 
that redundant parameters are grouped into clusters, de- 
noted by C1-C7, that to a large extent correspond to 
their topological localisation. Cluster CI contains pa- 
rameters describing receptor activation and signalling; 
C2: A20 synthesis and degradation; C3: IKK kinase 
post-translational modifications, InBa synthesis, degra- 
dation, phosphorylation, and interaction with IKKK and 
NF-/d3; C4: I^Ba synthesis, degradation and interaction 
with NF-ftB; (75: molecule numbers and cell character- 
isation (except ka20 and fc2); C6: NF-^B - DNA in- 
teractions; C7: nuclear shuttling. We can also consider 
two bigger clusters, one composed of elements of CI and 
C2, and the second of (C3-C8), which correspond to the 
external feedback loop controlled by A20, and internal 
feedback loop controlled by primary inhibitor IkBa, re- 
spectively. Parameters within the modules have a similar 
impact on model dynamics and their perturbations can 
be compensated by changes of other parameters within 
the module. Module by module compensation is less effi- 
cient i.e. correlation (mutual information) between clus- 
ters is smaller than between parameters within the clus- 
ters. 

As a second scenario we consider the behaviour of the 
system, X, to be defined entirely by the normalised nu- 
clear concentration of the NF-^B with measurement er- 
ror. The corresponding dendrogram is presented in Fig- 
ure 1A of SI. The redundancy between parameters is 
much stronger in this case (many low linkages) and the 
structure appears to be more random what reflects much 
lower information content of measurements. 
Robustness properties of the pathway. The dendrogram 
in Figure [4] gives a good characterisation of robustness 
of the response to TNF-a stimulation. In the first case 
(all variables define model behaviour) robustness is a 
topologically-local property i.e. system is robust to per- 
turbations that compensate each other locally. In the sec- 
ond scenario (nuclear NF-^B defines model behaviour), 
the system is robust to a much wider class of perturba- 
tions as most of the parameter pairs are mutually com- 
pensative. Below we also show that redundancy structure 
is stimulus/experiment specific. Stimuli that break func- 
tional redundancy are desired to infer model parameters 
from experimental data and we show that it is possible 
to design stimuli which reduce the redundancy. 
Redundancy and identifiability in available experimental 
data. Highly redundant parameters have almost identi- 
cal impact on observed experimental data and therefore 
render non-identifiability. There is a number of formal 
definitions of identifiability [TTJ [35j [36] . Parameters are 
structurally non-identifiable if the system exhibits iden- 
tical behaviour for two different values of these param- 
eters. Locally, this can be detected solely based on the 



FI. Structurally identifiable parameters can be practi- 
cally non-identifiable when for a given experimental data 
set likelihood is almost flat with respect to these pa- 
rameters [11 . In our setting this corresponds to the 
case where mutual information is high (alternatively p 
is close to 1) or the corresponding sensitivity coefficient 
is small. The existing definitions together with our ap- 
proach motivate to define identifiability in terms of re- 
dundancy. We introduce definition of (<S, £)-identifiably. 
We say that an element 6i of the vector = (0i,...,0fc) 
is (£, (^-identifiable if p{0i,6-i) < 1 — 5 (equivalently 
I{0i,0-i) < log(l/V<K2 -£))) and FI U (6) > ( (see SI), 
where 0_^ denotes all elements of the vector except 0$. 
This definition is rooted in the conventional statistics. 
If all other parameters were known, the standard devi- 
ation of a most efficient estimator of 0^, sd(0i), is given 
by the asymptotic formula sd(6i) = 1/ y/FIa(0). The 
( requirement demands the individual standard devia- 
tions to be smaller then l/vT- The S condition requires 
the asymp totic stan dard deviation not to increase more 
than (1/^/6(2 — 5))-fold when all elements of are es- 
timated at once compared to the previous case (see SI 
for detailed explanation). Using the constructed crite- 
ria we asked how many parameters of the model [T3l [30] 
can be estimated from data available in the literature. 
We select 9 papers that contain rich data sets on the 
dynamics of the NF-^B system [22-30 , which could be 
used for parameter inference. All experimental measure- 
ments are summarised in Table 1 in SI. We arbitrarily 
set S = 0.05, which corresponds to the increase of the 
asymptotic variance less than 10 times; and £ = 1. As we 
use logarithms parametrisation, i.e. log(0^) instead of 0^, 
setting £ = 1 corresponds to learning a parameter with 
an order of magnitude error, if all other parameters were 
known. Under this assumptions we found that 26 param- 
eters can be estimated. The identifiable parameters are 
plotted black in the Figure [5] Among non-identifiable 
parameters we found M, KN, KNN describing levels of 
receptors, IKKK kinase and IKK kinase respectively and 
parameters k2 and ka20 describing signalling mediated 
by receptors, IKKK and IKK. 

Informative future experiments. The redundancy struc- 
ture for all experiments of [22H30] presented in Fig. [5] 
indicates that certain parameters cannot be inferred as 
a result of their redundancy with other parameters of 
the model. Non-identifiabilty does not result from pa- 
rameters being insensitive individually. To find out if 
and to what extent the redundancy can be eliminated 
by designed stimuli we randomly searched a space of po- 
tential new stimuli (see SI) and selected 10 experiments 
that collectively give the highest number of identifiable 
parameters. To ensure practical value of our guidance 
we assumed that only the following entities were mea- 
sured: InBa protein (blotting), InBa mRNA, nuclear 
NF-ftB (fluorescence microscopy), and activity of IKK. 
The measurements were assumed to be taken every 5 
minutes for 160 minutes with normally distributed er- 
ror with standard deviation equal to the squared root of 
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FIG. 4. Functional redundancy in the NFftB signalling pathway. (A) Structure of the network. (B) Dendrogram of functional 
redundancy represented as dendrogram. Structure was computed assuming all model variables normalised by their maxima 
are observed with normal zero mean, unit standard deviation measurement error. Parameters were colour coded to represent 
biochemical elements of the network. Specific TNF-a stimulation to reflect physiological conditions (increase, plateau, decrease) 
was applied (see Fig. IB in SI). Clusters CI — C8 exhibit good correspondence with functional components of the network. 
A corresponding analysis under the assumption that the nuclear NF-^B trajectory is the only observable is presented in the 
Figure 1A in SI. (C) Normalised Fisher Information Matrix used to construct dendrogram. 



possible maximum of the measurement. We found that 
by caring out 10 experimental stimulations with TNF-a 
presented in Fig. [6] and Table 4 in the SI we can identify 
values of 7 new parameters. New experiments involve 
pulses with frequency ranging between 1/4 — 1/2 of the 
NF-ftB oscillation frequency which is approximately 100 
min. Such protocols have not been performed so far. 



III. DISCUSSION 

Tools to understand the relationship between param- 
eters and model dynamics seem to be of high relevance 
to utilise the potential of mathematical modelling within 
bimolecular sciences. The complexity of dynamic quanti- 
tative biology models makes their manipulation a subject 
of time consuming and laborious investigations. It is so 
because components and parameters are not independent 
(orthogonal) but exert their impact jointly through the 
network of interactions. A number of studies reported 



on the analysis of multi-parameter models. A notion of 
the sloppy model was introduced in [2j [15] [16] to de- 
scribe the model property of being sensitive only to a 
small number of linear combinations of parameters. New 
technique for sensitivity analysis that takes into account 
varying contributions of parameters into sensitivity spec- 
trum was proposed in [1 j. Information theory was ap- 
plied in [20] to reveal higher order interactions between 
model parameters. Novel methods proposed in [TTJ [21] 
can identify linear and non-linear relations between pa- 
rameters and detect non-identifiability in experimental 
data. A conjectured solution to the problem of over- 
parametrized models is to find and manipulate orthog- 
onal model components. This however creates abstract 
objects that may not have a relevant biological interpre- 
tation. Our method enables visualisation of similarities 
(redundancies) to find most redundant and most orthogo- 
nal biological components. It constitutes a unique mathe- 
matical framework to describe a variety of compensation 
like phenomena. Generality of the method enables its 
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FIG. 5. (A) Redundancy structure indicates identifiable parameters given experimental data published in [22-30 ]. Setting 
5 = 0.05, C — 1 we found that 26 parameters can be identified (plotted in black, the non-identifiable parameters in red. 
Sensitivity coefficients i.e diagonal elements of FI matrix are shown below the dendrogram. The heatmap plot presents the 
normalised FI matrix, its eigen spectrum is shown below. Linkages plotted in bold violate the condition p(6i, 0-i) < 1 — S. (B) 
Same as in (A) but for the experiments published in [22H30] together with 10 best experimental protocols (Fig. [6} found in a 
random search. 
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FIG. 6. TNF-a stimulation protocols designed to increase the 
number of identifiable parameters. Each row describes one of 
the 10 found temporal profiles The protocols together with 
data of [22H30] allow to estimate 7 new parameters compared 
to using data of 22-30 only. Details of the protocols are 
presented in the Table 4 in the SI . 



applicability to stochastic and deterministic models (see 
SI). A most detailed insight is provided if model param- 
eters are known. Nevertheless, information about redun- 
dancies can also be inferred directly from experimental 
data (see SI). The accuracy of inferred similarities de- 
pends on quality of available data. 

In the paper we focused on the computationally least 
demanding scenario of deterministic models with avail- 
able parameter estimates. We demonstrated the po- 
tential of our method by addressing relevant questions 
pertaining to the dynamics of the NF-/d3 system. We 
shown that functionally related parameters are topolog- 
ically co-localised. Modules formed by redundant pa- 



rameters impact the dynamics of the system in a more 
independent manner. This has consequences for model 
robustness and parameter identifiability, particularly im- 
plies that it is most difficult to infer parameters that are 
close in the networks topology. We also examined the 
literature available experimental protocols to show how 
redundancy disables inference of parameters of the NF- 
ftB dynamical model and suggested TNF-a stimulation 
protocols to break the redundancy and infer more pa- 
rameters. In the SI we shown how the method can be 
applied to analyse stochastic models of moderate size. 
Utilisation of the method for scenarios where only vague 
priors on parameter values are available requires further 
development of computational techniques. 
The introduced concept of redundancy shows how com- 
ponents of biological systems interrelate exerting the 
joint impact on observed biochemical dynamics. There- 
fore it has a tangible potential to overcome some of the 
difficulties resulting from the complexity of models in 
quantitative biology. 
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